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/  SUMMARY 

The  satellite  1971-10B  passed  through  exact  15th-order  resonance  on 
1981  March  30,  and  orbital  parameters  have  been  determined  at  52  epochs 
from  some  3500  observations  using  the  RAE  orbit  refinement  program,  PROP, 
between  1980  September  and  1981  October. 

The  variations  in  inclination  and  eccentricity  during  this  time  have 
been  analysed,  and  six  lumped  15th-' ^der  harmonic  coefficients  and  two 
30th-order  coefficients  have  been  evaluated.  The  15th-order  coefficients 
are  the  best  yet  obtained  for  an  orbital  inclination  near  65®;  and  ~r7 

previously  there  were  no  30th-order  coefficients  available  at  this 
inclination. 

The  lumped  coefficients  have  been  used  to  test  the  Goddard  Earth 
Model  GEM  10B:  there  is  good  agreement  for  seven  of  the  eight  coefficients 
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J  INTRODUCTION 

Cosmos  394  rocket,  designated  1971  —  1  OB ,  entered  orbit  on  1971  February  -9  and 
initially  had  the  following  orbital  elements1:  inclination  65.84°;  nodal  period  96.43  min; 
apogee  and  perigee  heights,  612  km  and  564  km  respectively;  and  eccentricity  0.003.  The 
satellite  will  decay  in  1984. 

By  1981  March  the  orbital  period  of  the  satellite  had  decreased  to  95  min  and  the 
orbit  was  passing  through  the  condition  of  15th-order  resonance.  This  occurs  when  the 
track  over  the  Earth  repeats  after  15  revolutions,  ie  when  the  satellite  makes  15  revolu¬ 
tions  while  the  Earth  spins  once  relative  to  the  orbital  plane.  For  this  Report  the 
orbit  of  the  satellite  has  been  determined  from  radar  and  visual  observations  between 

1980  September  and  1981  October,  using  the  RAE  orbit  refinement  program  PROP  in  the 

.  2 

PROP  6  version  .  Then  the  inclination  and  eccentricity  were  analysed  over  this  time  to 
evaluate  lumped  harmonic  coefficients  of  order  15  and  30  in  the  geopotential. 

The  analysis  of  this  orbit  at  15th-order  resonance  is  of  particular  interest  because 

of  its  inclination,  65. 8°.  In  previous  evaluations  of  the  individual  15th-order 

.3 

harmonic  coefficients  in  the  geopotential  from  the  lumped  harmonic  coefficients  derived 
from  the  analysis  of  satellite  orbits  at  various  inclinations,  there  have  been  no 
accurate  lumped  coefficients  available  for  an  inclination  near  65°. 

2  THE  OBSERVATIONS,  ORBITS  AND  OBSERVATIONAL  ACCURACY 

2.1  Observations 

The  orbit  of  197 1-1 0B  has  been  determined  from  some  3500  observations,  between  1980 
September  3  and  1981  October  28,  at  52  epochs.  The  observations  used  fell  into  three 
groups  and  a  breakdown  of  the  number  and  type  used  on  each  orbit  is  given  in  Table  I . 

The  first  group  consists  of  visual  observations  made  by  volunteer  observers  who 
report  to  the  University  of  Aston  in  Birmingham  and  these  observations  usually  have 
accuracies  between  1  and  4  minutes  of  arc.  The  second  group  is  made  up  of  Navspasur 
observations  kindly  supplied  by  the  US  Naval  Research  Laboratory,  with  accuracies  of 
about  2  minutes  of  arc.  The  third  group  comprises  kinetheodolite  observations  made  at 
the  South  African  Astronomical  Observatory,  accurate  to  about  1  minute  of  arc.  Observa¬ 
tions  from  the  kinetheodolite  were  only  available  on  orbits  1-18,  as  the  instrument  closed 
down  at  the  end  of  January  1981. 

2.2  The  orbits 

The  orbits  were  determined  at  approximately  8-day  intervals  from  1980  September  to 

1981  October  with  the  aid  of  the  RAE  orbit  refinement  program  in  the  PROP  6  version,  and 
the  orbital  elements  at  each  of  the  52  epochs  are  listed  in  Table  2,  with  the  standard 
deviations  below  each  value.  The  epoch  for  each  orbit  is  at  00  hours  on  the  day  indicated, 
and  the  PROP  program  fits  the  mean  anomaly  M  by  a  polynomial  of  the  form 

M  -  M.  +  M.t  ♦  M,t2  +  M,t3  ♦  M.t4  +  M.t5  ,  (1) 

U  12  3  4  :>  * 

where  t  is  the  time  measured  fiom  epoch,  and  the  number  of  M  coefficients  used  depends 


Table  1 

Sources  of  observations  used  in  each  orbit 


Orbit 

Ho. 

Source  of  observation 

Orbit 

Source  of  observation 

Visual 

US 

Navy 

Cape 

kine 

Total 

No. 

Visual 

US 

Navy 

Total 

1 

51 

51 

27 

17 

76 

93 

2 

8 

78 

66 

28 

1 

75 

76 

3 

7 

61 

68 

29 

4 

59 

63 

4 

9 

74 

83 

30 

55 

55 

5 

8 

59 

67 

31 

1 

46 

47 

6 

7 

61 

68 

32 

4 

41 

45 

7 

6 

73 

79 

33 

7 

71 

78 

8 

45 

3 

48 

34 

4 

88 

92 

9 

8 

66 

74 

35 

65 

65 

10 

1 

42 

1 

44 

36 

5 

59 

64 

11 

3 

69 

2 

75 

37 

6 

50 

56 

12 

4 

84 

88 

38 

69 

69 

13 

49 

49 

39 

56 

56 

14 

16 

54 

70 

40 

4 

62 

66 

15 

8 

55 

63 

41 

62 

62 

16 

18 

80 

98 

42 

54 

54 

17 

4 

64 

68 

43 

11 

53 

64 

18 

1 

55 

9 

65 

44 

18 

67 

85 

19 

3 

47 

50 

45 

57 

35 

92 

20 

6 

48 

54 

46 

4 

53 

57 

21 

48 

48 

47 

4 

50 

54 

22 

86 

86 

48 

15 

56 

71 

23 

2 

113 

115 

49 

7 

45 

52 

24 

6 

68 

74 

50 

4 

63 

67 

25 

3 

75 

78 

51 

78 

78 

26 

10 

80 

90 

52 

3 

44 

47 

*  on  the  drag.  For  197 1-1  OB,  which  was  in  a  nearly  circular  orbit  at  about  500  km  height, 

•  Mq  -  Mj  were  the  coefficients  required  on  30  orbits;  but  for  the  remaining  22  orbits 
only  the  coefficients  Mq  -  Mj  were  needed. 

The  value  of  c  ,  the  parameter  indicating  the  measure  of  fit  of  the  observations 
to  the  orbit,  varied  between  0.33  and  1.16,  and  had  an  average  value  of  0.62. 

_  For  all  52  orbits  the  standard  deviations  in  inclination  were  between  0.0007°  and 

0.0024°,  the  average  being  0.0013°  which  is  equivalent  to  about  150  m  in  distance.  The 
values  of  eccentricity  had  an  average  standard  deviation  of  0.000013  over  the  52  orbits, 
equivalent  to  90  m  in  distance.  The  individual  standard  deviations  in  e  varied  between 
0.000006  and  0.000027.  The  right  ascension  of  the  node  had  nearly  the  same  accuracy  as 
the  inclination. 

V  ♦' 

2.3  The  accuracy  of  the  observations 

The  residuals  of  the  observations  have  been  obtained  using  the  ORES  computer 
program  and  sent  to  the  observers.  Table  3  lists  the  accuracies  of  the  observing 
stations  with  five  or  more  observations  accepted  in  the  orbit  determinations.  The  Table 
shows  the  number  of  observations  used  from  a  particular  station,  divided  into  those 
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Values  of  the  orbital  parameters  at  52  epochs  with  standard  deviations 
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accepted  and  rejected  in  the  orbit  determinations.  The  US  Navy  observations  from 
station  29  are  geocentric,  and  if  they  were  given  in  the  same  form  as  the  topocentric  • 

observations,  their  angular  rms  residuals  would  increase  by  a  factor  of  about  five. 

The  residuals  for  station  29  were  similar  to  their  usual  standard.  Those  for 
stations  1-6  were  slightly  larger,  but  this  is  because  low-elevation  observations  (<20°) 
were  retained  with  their  assumed  error  doubled.  The  kinetheodolite  observations  had  an 
average  accuracy  of  1.8  minutes  of  arc  and  the  visual  observations  averaged  about 
5  minutes  of  arc.  These  accuracies  are  slightly  worse  than  usual,  for  no  obvious  reason. 

Table  3 

Residuals  for  observing  stations  with  five  or  more  • 

observations  accepted  in  the  final 
orbit  determination 


Station 

Number  of 
observations 

Rms  residuals 

Accepted 

Rej ected 

Range 

km 

Minutes  of 

arc 

RA 

Dec 

Total 

1  US  Navy 

315 

47 

2.6 

2.9 

3.9 

2  US  Navy 

216 

32 

3.3 

3.0 

4.5 

3  US  Navy 

166 

23 

4.3 

3.7 

5.7 

4  US  Navy 

186 

33 

4.0 

3.2 

5.1 

5  US  Navy 

228 

21 

2.6 

2.8 

3.3 

6  US  Navy 

266 

19 

2.7 

2.1 

3.4 

29  US  Navy 

1314 

35 

0.6 

0.3 

0.4 

2183  Pengwern 

6 

1 

12.7 

9.2 

15.7 

2265  Farnham 

14 

3 

4.5 

6.3 

7.7 

2392  Cowbeech 

5 

0 

1.4 

2.2 

2.6 

2414  Bournemouth 

56 

0 

3.2 

3.4 

4.7 

2418  Sunningdale 

23  • 

1 

4.1 

4.1 

5.8 

2420  Willowbrae 

85 

29 

3.5 

2.6 

4.4 

2577  Cape  kinetheodolite 

12 

3 

1.3 

1.2 

1.8 

4156  Apeldoorn 

9 

1 

2.9 

2.8 

4.1 

3  THE  EQUATIONS  AT  1 5TH-0RDER  RESONANCE 

The  theory  for  general  8:a  resonance  has  been  given  before,  in  Ref  3,  where  all 
parameters  used  in  the  following  equations  are  defined. 

The  theoretical  equations  at  15th-order  resonance  for  the  variation  of  inclination 
and  eccentricity  are  as  follows: 


84008 
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di  „  n  /R  V 5 1 
dc  sin  i  \a  /  I 


0,1 


(15  -  cos  i>F|5fI5f7jC15  6in^-SJ5  cos  4> 

J,  o 


T*  (15)(l)F.6,.5,8jS15  sin(*  -  “>  +  Cl 5  C°8($  '  u)| 

T*  (15  -  2  cos  i)(|)F,6fl5>7js,5’  sin(*  ♦  *>  ♦  CJ5’  cos(*  ♦  u)J 


+  terms  m 


(|£e)  " 1  cos 


(k!) 


.  (yt  -  qu)> 


sin 


(2) 


and 


de 

dt 


-  17 


♦  13 


8  (S^  sin(*  -  in)  *  Cjj  cos($  -  o>> 


n/Rl'5  -  /p0*1  .  4  A 

l(l)  <*15,15,7^.5  Sln*'S15  cos 

/R\-  t_*»0 

7(l)F16,15, 

»(l)Vl5.7  j'  1,2 

e- -  jq  -  l(k  +  q)e2}  “Jw  -  qui) 


J,0 

:15 

--*.2 


{S)5  sin(«  +  u)  *  CJ5  cos(0  +  m) 


+  terms  in 


(lqi)! 


The  resonance  angle  f  is  given  by 


w  *  M  +  1 5(n  -  v)  , 


(3) 


(4) 


where  v  is  the  sidereal  angle,  and  at  exact  resonance  $  «  0  .  In  equations  (2)  and  (31 
only  the  three  terms  with  (Y,q)  B  < 1 , 0) (1,1)  and  (1,-1)  are  given  explicitly. 

The  three  pairs  of  lumped  coefficients  in  equations  (2)  and  (3)  can  be  expressed 
as  linear  sums  of  individual  15th-order  coefficients 


.q*k 


Urn 


_q»k 

and  S 

m 


£<!?• 


£m 


(5) 


and  these  are  given  explicitly,  with  the  expressed  in  terms  of  the  F  inclination 

functions,  ir.  Ref  3. 

» 

The  orders  of  magnitude  of  the  (y,q)  terms  in  the  equation  for  di/dt  can  be 
•  —  _  „c  o 

estimated,  as  the  (or  S  )  are  expected  to  be  of  order  10  /l  ,  so  the  value 

-q,k  _q,k  , 

o*  Cm  (°r  Sm  ^  can  be  taken  to  °f  order  *  10  ’  /  •  The  F  and  Q 

coefficients  can  be  calculated  for  specified  values  of  (y»q)  using  the  computer  program 


PROF.  The  expected  orders  of  magnitude  are  given  in  Table  4  without  the  factor 


sin 


(!f 


(15  -  cos  i)  * 


Table  4 


Estimated  orders  of  magnitude  of  terms  in 
equation  for  di/dt 


Y 

<1 

1 

2 

3 

0 

19 

2.4 

0.4 

+1 

0.4 

0.1 

- 

-1 

0.3 

0.1 

The  expected  orders  o£  magnitude  for  the  (y.q)  terms  in  the  equation  for  de/dt 

can  be  estimated  in  a  similar  manner  to  di/dt  and  they  are  given  in  Table  5  without 
/r\J5  —9 

the  factor  n(— j  *  10  . 


Table  5 


I 


Estimated  orders  of  magnitude  of  terms  in 
equation  for  de/dt 


1 

2 

3 

0 

0.02 

- 

- 

+  i 

157 

21 

4 

-i 

136 

18 

3 

+2 

3 

- 

- 

-2 

2 

The  orders  of  magnitude  of  the  (y,q)  terms  in  Tables  4  and  5  suggest  that  taking 
the  (y,q)  "  0,0)  (2,0)  terms  for  inclination  will  only  introduce  an  error  of  27.  or  less, 
and  taking  the  (y,q)  “  (1,±1)(2,±1)  terms  for  eccentricity  only  3Z  or  less. 

4  THE  ANALYSIS  OF  INCLINATION  AT  I 5TH -ORDER  RESONANCE 

The  satellite  1971  —  1  OB  passed  through  exact  15th-order  resonance  on  1981  March  30. 
The  variation  of  the  resonance  angle  $  ,  given  by  equation  (4),  and  the  rate  of  change 
of  the  resonance  angle,  4  ,  whiefti  increased  from  -18  to  +18  deg/day  between  1980 
September  and  1981  October,  are  shown  in  Fig  1. 

Before  the  changes  in  inclination  due  to  resonance  can  be  analysed,  all  other  known 
perturbations  must  be  removed.  The  52  values  of  inclination  available  in  Table  2  were 
cleared  of  lunisolar  and  zonal  harmonic  perturbations,  by  using  the  computer  program 
PROD^  with  l-day  integration  steps,  and  the  perturbations  due  to  the  J 2  ,  tesseral 
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harmonic,  which  is  recorded  on  each  PROP  run,  were  also  removed.  Perturbations  due  to 
earth  and  ocean  tides^  should  not  exceed  50  m  and  need  not  be  considered  here  as  the 
values  of  inclination  have  average  standard  deviations  of  150  m. 

The  values  of  inclination,  cleared  of  the  above  perturbations,  and  with  the 
standard  deviations  quoted  in  Table  2,  were  fitted  with  equation  (2)  in  integrated  form 
using  the  THROE^  computer  program  which  removes  perturbations  due  to  atmospheric  rotation 
and  lunisolar  precession  of  the  Earth's  axis.  The  density  scale  height  H  was  taker,  as 

g 

73  km,  appropriate  to  a  height  of  512  km,  0.2H  above  perigee  ,  and  the  atmospheric  rota- 

9 

tion  rate,  A  ,  was  taken  as  0.8  rev/day.  The  values  of  M2  were  altered  to  mean 
values,  M2  •  by  the  technique  described  in  Ref  10. 

The  first  THROE  fitting  of  i  ,  with  the  (y,q)  =  (1,0)  terms  only,  gave  a  value  of 
e  ,  the  parameter  indicating  the  measure  of  fit,  of  0.726  and  the  values  of  the  lumped 
l5th-order  harmonics  were: 

q  °»>  q-0.1 

10  C]5  *  2.6  i  3.9  10  s15  =  3.1  i  °-9  •  (6) 

—5  2 

The  values  of  these  coefficients  were  surprisingly  small:  the  10  /l  rule  indicated  that 

-9 

they  should  be  of  order  81  *  10  .  This  means  that  the  (y,q)  «  (1,0)  terms  in  Table  4 

are  only  about  1/30  of  the  expected  value:  it  was,  therefore,  necessary  to  include  the 

_ 0  2 

(2,0)  terms  in  the  fitting.  The  (C,S),I  coefficients  in  the  (2,0)  terms  should  be 

“9  2  0  1 

of  order  52  *  10  according  to  the  10""Vj (  rule  and  as  the  (C,S)jj  terms  are  so 

small,  the  (2,0)  terms  may  be  dominant. 

The  second  fitting,  with  (y,q)  =  (1,0) (2,0)  terms  included,  gave: 

q  q_°.> 

'°C15  -  -  0.7  r  4.1  ,09S)5  =  2.4:  3.9 

0  2  0  2 

q_  q_  * 

10  C  =  -  54  :  27  10  S30  =  59  -  40 

.0,2  _0,2 

with  t  =  0.705.  Since  the  values  of  C,,.  ,  art  slightlv  greater  than  the  10  '  !'■ 

-0,1  .0,1  30  30 

rule  and  the  ,  S  j  are  only  0.03  times  the  expected  value,  the  relative  magnitudes 

of  the  (1,0)  and  (2,0)  terms  in  equation  (7)  are  not  19  and  2,4  respectively,  as  in 
Table  4,  but  0.6  and  2.7.  Thus  the  (2,0)  term  is  the  more  important  and  the  solution  for 
(1,0)  alone,  equation  (6),  cannot  be  accepted. 

The  values  of  inclination  were  next  fitted  with  (y,q)  =  (1,0) (2,0)  and  (3, 

This  fitting  was  not  so  good:  the  (3,0)  terms  were  undetermined  and  the  standard  devia¬ 
tions  on  the  other  four  coefficients  were  increased.  Omitting  the  small  (1,0)  terms  did 
not  improve  the  (3,0)  terms,  as  they  were  still  undetermined.  A  fitting  with  just 
(y,q)  «  (2,0),  whilst^g jving  slightly  lower  standard  deviations,  was  not  considered 
appropriate  as  the  coefficients,  although  small,  have  standard  deviations  c: 

the  same  magnitude  as  1 5th-order' coefficients  at  other  inclinations^  and  therefore 
provide  useful  information  on  15th-order  harmonics.  So  the  values  in  equation  (“■  vert 
preferred. 
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The  effects  of  changing  the  values  of  H  and  A  were  investigated.  The  altera¬ 
tion  of  H  had  no  effect  on  the  value  of  e  or  the  lumped  coefficients.  Changing 
the  value  of  A  to  0.7  or  0.9  increased  the  value  of  e  ,  thus  confirming  the  choice 
of  the  value  0.8.  The  effect  of  shortening  the  run  was  also  investigated,  by  dropping 
the  first  five  and  the  last  five  values  of  inclination.  This  resulted  in  a  value  of  e 
only  slightly  lower,  and  the  values  of  the  harmonic  coefficients  were  not  significantly 
altered.  As  the  change  in  the  inclination  at  15th-order  resonance  is  very  small,  it 
seemed  preferable  to  keep  as  many  points  as  possible  on  either  side  of  exact  resonance 
to  firmly  establish  the  values  of  inclination  before  and  after  resonance.  So  the 
coefficients  obtained  using  all  52  values  of  inclination  were  preferred. 

The  values  of  inclination,  cleared  of  all  known  perturbations  except  those  due 
to  resonance,  are  plotted  in  Fig  2,  with  the  theoretical  curve  derived  from  the  THROE 
fitting  with  the  (y,q)  “  (1,0)(2,0)  terms,  which  gave  the  coefficients  in  equation  (7). 

Although  the  change  in  inclination  as  a  result  of  passing  through  resonance  is  very 
small,  0.0010°,  Fig  2  shows  that  the  fitting  of  the  values  is  completely  satisfactory. 

The  values  of  the  15th-order  lumped  coefficients  in  equation  (7)  may  look  inaccurate,  but 
their  standard  deviations  correspond  to  errors  in  geoid  height  of  only  about  2  cm. 

5  THE  ANALYSIS  OF  ECCENTRICITY  AT  1 5TH-0RDER  RESONANCE 

The  perturbations  in  eccentricity  due  to  zonal  harmonics  and  atmospheric  drag,  and 
the  (much  smaller)  lunisolar  perturbations,  all  have  to  be  removed  before  the  changes  due 
to  resonance  can  be  analysed. 

The  removal  of  the  zonal  harmonic  and  lunisolar  perturbations  was  performed  by 
using  the  PROD  computer  program.  The  variation  of  e  with  time  on  a  drag-free  orbit 
was  calculated  by  numerical  integration  at  )-day  intervals  with  the  aid  of  this  program, 
and  this  drag-free  orbit  gives  the  changes  in  e  due  to  zonal  harmonic  and  lunisolar 
perturbations.  In  Fig  3a  the  values  of  eccentricity  from  Table  2  are  plotted  as  circles 
and  the  triangles  show  the  values  of  e  after  removal  of  zonal  harmonic  and  lunisolar 
perturbations,  e  -  AeZHL  •  The  remaining  variation,  indicated  by  the  triangles,  is  due 
to  air  drag  and  resonance. 

The  air  drag  model  within  THROE  takes  no  account  of  the  day-to-night  variation  in 

density,  which  is  important  for  1971-10B.  So  the  correction  due  to  drag  in  an  atmosphere 

with  day-to-night  variation  was  calculated  using  the  theory  given  in  the  Appendix  of 

Ref  11,  with  scale  height  H  taken  as  74  km,  appropriate  to  a  mean  height  of  521  km, 
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which  is  0.2H  above  perigee  .  This  theory  assumes  that  the  atmospheric  density  depends 

on  the  geocentric  angular  distance  from  the  point  of  maximum  density,  taken  to  be  on  the 

1  3 

equator  at  either  14  h  local  time  (the  time  indicated  by  CIRA  197Z)  or  at  16  h  local 

;  j  4 

time  (the  time  indicated  in  Jacchia’s  1977  model  ).  The  correction  Ae^  for  drag  with 
day  -to-night  variation  is  shown  in  Fig  3b  for  local  times  of  14  h  and  16  h. 

The  52  values  of  eccentricity  from  Table  2  could  then  be  cleared  of  perturbations 
due  to  zonal  harmonics,  lunisolar  forces  and  atmospheric  drag,  and  the  resulting  values 
of  e  should  show  a  variation  due  to  resonance  alone.  These  values  of  e  were  then 


fitted  with  equation  (3)  in  integrated  form,  using  THROE.  The  program  was  run  with  no 
adjustment  being  made  to  e  to  compensate  for  zonal  harmonic  and  atmospheric  perturba¬ 
tions  as  the  values  of  e  had  already  been  cleared  of  these  perturbations. 

The  first  THROE  run  was  with  the  values  of  e  changed  on  the  assumption  that  the 
point  of  maximum  density  was  at  14  h  local  time  and  the  second  run  at  16  h  local  time. 

The  first  run  gave  e  «  3.61  and  the  second  e  ■  2.98;  so  the  values  of  e  adjusted  for 
16  h  local  time  were  adopted.  The  last  six  values  of  e  were  ill-fitting,  however, 
and  when  they  were  omitted,  a  further  fit  with  THROE  gave  a  value  of  t  reduced  by  7 Z 
to  2.76.  These  fittings  were  with  (y,q)  ■  (l,l)(l,-|). 

Previously1^,  when  fitting  values  of  eccentricity  with  THROE  at  inclinations  near 
the  critical  inclination  of  63.4°  it  was  found  that  some  adjustment  of  the  zonal  harmonics 
was  needed.  This  arose  because,  when  the  set  of  odd  zonal  harmonics  used  in  THROE  were 
evaluated16,  there  were  only  a  few  satellites  available  for  analysis  with  inclinations  in 
this  critical  region.  Therefore  the  values  of  the  odd  zonal  harmonics  used  in  THROE  nay 
not  be  adequately  modelling  the  variation  in  eccentricity  for  1 971  —  1  OB ,  which  is  at  an 

inclination  of  65.83°.  The  alterations  were  made  by  adding  an  increment  AJ.,  to  J.  , 

“6  ^  ** 
which  has  a  value  of  -2.531  *  10  in  the  THROE  model.  Fig  4  shows  the  variation  of  t 

with  AJ3  .  The  lowest  value  of  e  ,  2.751,  was  attained  with  AJ^  -0.013  ,  and  this 

was  accepted  as  the  best  fitting.  This  implies  that  the  lumped  harmonic  Y  defined  in 

Ref  17  is  in  error  by  0.013:  since  the  errors  in  the  values  of  Y  implicit  in  THROE16  at 

nearby  inclinations  are  all  at  least  0.03,  this  refinement  is  well  within  the  limits  of 

error  of  the  model. 

This  final  fitting  of  the  values  of  eccentricity  with  THROE  (with  AJ^  “  -0.013) 
yielded  values  of  the  lumped  harmonics  as  follows: 

q_*  *0  q_* 

10  C)5  *  35.1  -  11.7  10  s15  *  -  !3*8  11  n-° 

-1,2  -1  2 

109c)5  -  -20.0  c  10.7  109S,5’  *  -  18.9  ±  10.3 

The  values  of  eccentricity,  cleared  of  all  known  perturbations  except  those  due  to 
resonance,  are  plotted  as  circles  in  Fig  5  with  the  standard  deviations  from  Table  2, 

The  theoretical  curve  is  that  of  the  final  THROE  run,  which  gave  the  coefficients  in 
equation  (8).  A  fitting  of  e  by  THROE  with  (y,q)  m  (l,il)  and  (2,±1)  terms  was  also 
tried,  but  this  resulted  in  a  larger  value  of  c  and  undetermined  coefficients. 

The  fitting  of  eccentricity  in  Fig  5  is  far  from  satisfactory,  as  the  high  value 

of  e  (2.75)  indicates.  Between  the  initial  point  at  MJD  44485  and  MJD  44625  there 

appear  to  be  lj  cycles  of  a  sinusoidal  oscillation  in  the  observational  values,  with 
period  about  85  days  and  amplitude  0.00008  ( ie  about  500  m) .  Fig  3b  shows  that  the  day- 
to-night  variation  in  drag  has  a  similar  periodicity,  but  its  effects  are  too  small  to 
cancel  completely  the  large  oscillation  visible  in  the  observational  values,  which  has  an 
amplitude  of  about  0.00012  (  ic  800  m) .  The  explanation  of  this  anomaly  remains  uncertain, 
but  it  seems  likely  to  be  attributable  to  inadequate  modelling  of  the  day-to-night 
variation  in  air  density,  perhaps  due  to  asymnetry  about  the  equator  at  times  near  the 
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solstices,  or  perhaps  resulting  from  a  real  atmospheric  anomaly  between  1980  September 
and  1981  January. 

Whatever  the  cause  of  the  anomaly  in  e  ,  it  is  fortunate  that  the  fitted  curve 
cuts  through  the  mean  of  the  oscillation,  and  should,  therefore,  still  provide  reliable 
values  of  the  coefficients,  even  if  the  standard  deviations  are  disappointingly  large. 

6  INCLINATION  AND  ECCENTRICITY  FITTED  SIMULTANEOUSLY 

The  values  of  inclination  and  eccentricity  fitted  separately  by  THROE  can  be 

|  fi 

fitted  simultaneously  using  the  computer  program  SIMRES  .  This  program  combines  the 
results  from  a  number  of  THROE  runs  (each  with  the  same  set  of  (y,q)  terms),  and  produces 
a  single  set  of  coefficients  to  fit  the  data.  For  this  SIMRES  fitting  the  results  of 
THROE  runs  with  (y,q)  *  ( 1 ,0) ( 1 , 1 ) ( 1 ,-l )  and  (2,0)  were  used.  The  SIMRES  program  allows 
a  choice  of  weighting,  so  that  the  contributing  THROE  runs  can  be  given  more  or  less 
weight  according  to  their  accuracy  of  fit,  which  is  indicated  by  the  value  of  c  . 

The  THROE  fitting  of  inclination  with  (y.q)  •  (1 ,0) ( 1 , 1 )( I ,-I )  and  (2,0)  gave 
e  «  0.602,  and  for  eccentricity  the  value  of  c  was  2.843,  when  fitted  with  the  same 
terms.  For  the  SIMRES  fitting,  therefore,  the  weighting  of  e  was  down-graded  by  a 
factor  equal  to  the  ratio  of  the  values  of  c  on  the  THROE  fittings,  namely  4.723 
(-2.843/0.602).  The  values  of  the  lumped  coefficients  given  by  this  SIMRES  fitting  are: 
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The  fittings,  if  shown  pictorially  in  Figs  2  and  5,  would  be  indistinguishable  from  the 
curves  given  by  the  THROE  fittings. 

The  lumped  coefficients  obtained  from  the  fitting  of  i  and  e  separately,  given 
in  equations  (7)  and  (8)  are  very  similar  to  those  obtained  from  the  simultaneous 
fitting  given  in  equation  (9).  So  either  set  of  values  could  be  used  in  future  deter¬ 
minations  of  the  individual  ISth-order  harmonic  coefficients  in  the  geopotential. 

7  THE  EQUATIONS  FOR  THE  INDIVIDUAL  COEFFICIENTS 

The  lumped  coefficients  can  be  expressed  as  linear  sums  of  individual  l5th-order 
coefficients,  see  equation  (5),  and  may  be  written 


-0.' 

_ 

0. 1- 

0.1- 

C,5 

Cl 5, 1 5 

*  QI7  C 1 7, 1 5  * 

Q1 9  C19, 1 5  *  *** 

(10) 

J.O 

CI5 

^16,15 

-  n' .0 7  + 

QI8  C18, 15 

Q2o  So, is  ■*■••• 

(11) 

ROOW 


15 


“1.2 

o”''2C 

QI8  18,15 

-1,2 

<15 

’  <16,15  4 

4«!8  ' 

_0,2 

-0.2= 

-0,2= 

<30 

*  <30,30  4 

Q32  <32,30 

4  <>34  < 

♦  ... 


♦  ... 


(12) 

(13) 


and  similarly  for  S  ,  on  replacing  C  by  S  throughout. 

The  Q  coefficients  in  equations  (10)  to  (13)  can  be  evaluated  using  the  computer 
program  PROF,  and  numerical  versions  of  equations  (10)  to  (13)  are  given  below.  In 
equations  (10)  to  (12)  the  Q  terms  included  go  as  far  as  Q^'.  Q^g0  and  Q3'*‘ 
respectively,  to  correspond  with  Ref  3,  and  in  equation  (13)  as  far  as  Q<,~  »  to 
correspond  with  Ref  19. 

‘“s'  m  <15,15  -  !-83°517,15  -  °‘35°<19.15  *  °-74l<2»,15  *  °*701<23,15  +  °'107<25,15 
-  °*348<27,15  -  °*387<29,15  ’  °*l40<31,15  4  °*I2,<33,15  +  °*215<35,15 
+  0-,36S7,15‘  0*007539,»5  04) 
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and  similarly  for  S  ,  on  replacing  C  by  S  throughout. 

8  LUMPED  COEFFICIENTS  FROM  1 97 1 - 1  OB  COMPARED  WITH  THOSE  OBTAINED  FRO*  GEM  10B 

The  values  of  individual .coefficients  from  the  comprehensive  geopotencial  model, 
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Goddard  Earth  Model  GEM  10B  ,  vhich  extends  to  degree  and  order  36,  may  be  substituted 

into  equations  (I4)-(17)  and  the  corresponding  S  equations  to  give  lumped  coefficients. 


Table  6 


Values  of  lumped  coefficients 
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The  15th-  and  30th-order  lumped  coefficients  thus  obtained  from  GEM  10B  are  given  in 
Table  6,  together  with  those  from  the  THROE  fittings  of  inclination  and  eccentricity  for 

3 

1 97 1  —  1  OB.  In  Table  6  the  accuracy  of  the  GEM  10B  individual  coefficients  is  assumed  to 

— q 

be  3  *  10  ,  and  the  accuracies  for  the  lumped  coefficients  have  been  calculated  on  this 

assumption. 

GEM  JOB  emerges  well  from  this  test  of  its  accuracy  by  a  completely  independent 
method.  The  two  sets  of  values  agree  well:  seven  of  the  eight  coefficients  differ  by 
less  than  the  sum  of  their  standard  deviations;  for  the  eighth  coefficient  (Cj*  )  the 
values  differ  by  twice  the  sum  of  their  standard  deviations. 

It  is  not  fair  to  compare  the  values  obtained  from  197 1-1  OB  with  those  from  other 
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comprehensive  models,  such  as  GRIM  3  or  the  1981  model  by  R.H.  Rapp  at  Ohio  State 

University,  because  both  these  models  incorporated  !5th-order  terms  from  Ref  3  in  their 

solutions.  The  values  of  lumped  coefficients  derived  from  satellites  at  inclinations 

near  65°  which  were  used  in  Ref  3  were  known  to  be  of  poor  accuracy,  so  any  comparison 

with  the  two  models  mentioned  would  only  show  where  the  solution  of  Ref  3  was  in  error. 

A  more  fruitful  procedure  is  to  use  these  values  from  1971—1  OB  for  revising  Ref  3;  and 
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this  has  now  been  done  , 

9  CONCLUSIONS 

The  orbit  of  1971—1  OB  has  been  determined  at  52  epochs,  from  some  3500  observations, 
during  the  time  when  the  orbit  was  affected  by  !5th-order  resonance.  The  average 
accuracy  of  inclination  and  eccentricity  on  the  52  orbits  was  equivalent  to  150  m  and 
90  m  in  distance  respectively. 

The  variations  in  inclination  and  eccentricity  have  been  analysed:  six  !5th-order 

and  two  30th-order  lumped  harmonic  coefficients  have  been  evaluated.  The  six  !5th-order 

lumped  coefficients  have  already  been  used  in  new  determinations  of  the  individual  15th- 

order  harmonic  coefficients  in  the  geopotential  from  lumped  harmonic  coefficients  derived 

23 

from  the  analysis  of  satellite  orbits  at  various  inclinations  .  Previously  there  were 
no  accurate  lumped  coefficients  available  at  an  inclination  near  65°  and  the  addition  of 
the  results  from  1 97 1  — 1 0B  has  greatiy  improved  the  solutions. 

a 
c 
c 

oc 


Though  the  values  obtained  here  for  1 97 1—1 0B  are  better  than  any  previous  results 
for  inclination  near  65°,  it  is  hoped  to  determine  still  better  values,  especially  for 
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even-degree  coefficients,  from  the  satellite  1965-14A  which  has  recently  been  passing 
through  15th-order  resonance  at  a  slower  rate. 

The  two  30th-order  lumped  coefficients  should  be  very  valuable  in  a  determination 

of  the  individual  30th-order  coeff icients,  because  previously  there  were  no  results 

o  23 

available  near  63  inclination.  This  expectation  has  already  been  confirmed  , 

The  eight  lumped  coefficients  obtained  from  1 97 1  —  I  OB  have  been  used  to  test  GEM  10B: 
good  agreement  was  obtained. 
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